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MATHEMATICAL MODEL OF WATER TRANSPORT IN BACON AND ALKALINE 
MATRIX-TYPE HYDROGEN-OXYGEN FUEL CELLS 
by Paul R. Prokopius and Robert W„ Easter 
Lewis Research Center 

SUMMARY 


As part of a continuing program in the study of the water transport processes in 
hydrogen- oxygen fuel cells, a mathematical model was developed that simulates the 
transport in both Bacon and alkaline - matrix cells. 

Aside from identifying and defining the processes controlling water transport in 
these two types of cells, the model appears to be sufficiently general to provide a basis 
for simulating the transport processes in other state-of-the-art fuel cells. 

The model specifically identifies the basic water transport process in matrix cells 
to be higher than first order. This is significant in that the electrolyte drying which 
occurs following a matrix cell load increase can be partially attributed to this higher 
order characteristic and thus is not completely the result of operating temperature 
changes as is commonly believed. 

The model was developed by applying the general laws of mass continuity and dif- 
fusive mass transfer to the transport zones in a lumped parameter mass transport anal- 
ysis. The result is a simulation that comprises five simultaneous, nonlinear differen- 
tial equations that predict the movement of both water and hydroxide ions within a cell. 

Validation of the model was accomplished by comparing analytically obtained water 
transport data with corresponding data obtained in previously conducted experimental 
water transport studies. 


INTRODUCTION 


Control of the water inventory in the electrolyte of hydrogen- oxygen fuel cells is 
the prime requisite in maintaining cell electrolyte concentration within design limits 
and hence contributes to safe and efficient operation during both normal and emergency 
operating periods. The design of systems needed to achieve this type of control could 


be greatly aided by an understanding of the mechanisms that regulate the movement and 
hence the rejection of product water from an operating cell. On this basis a program is 
being conducted at the Lewis Research Center to identify and mathematically simulate 
the processes controlling fuel- cell water transport for both steady- state and transient 
operation. 

The water removal technique studied in this program is one in which the product 
water is evaporated into a circulating reactant (Hg) stream. The rate of evaporation is 
controlled by regulating the parameters, particularly the inlet humidity, of the circulat- 
ing hydrogen stream. 

As a first step in gaining an understanding of fuel- cell water transport, experimen- 
tal and analytical studies of the water transport in a Bacon cell were conducted. The 
results of these studies are presented in reference 1. On completion of the Bacon cell 
work, the program emphasis was shifted to matrix cells, and an experimental dynamics 
study was conducted and reported on in reference 2. From the data obtained in this 
study, it became evident that the mathematical Bacon cell model lacked the generality 
to be applicable to a matrix cell. Therefore, using the theory developed in the 
Bacon cell model as a starting point, a model has been developed that is sufficiently 
general to simulate the water transport of both Bacon and alkaline matrix cells. 

As in the case of the Bacon cell, this model was developed on a lumped parameter 
basis. The movement of water within the cell is simulated by applying the equation of 
continuity and the laws of diffusive transport to the various zones or lumps into which 
the cell structure was schematically divided. 

The model equations were simulated on an analog computer on which the perturba- 
tions applied in the experimental water transport studies were reproduced analytically. 
By comparing the analytically obtained response data with the corresponding experimen- 
tal data, the model was found to be a valid simulation of both the Bacon and matrix cell 
water transport. In addition, the model appears to be sufficiently general to serve as 
the nucleus for the transport simulation for cells that employ methods other than the 
circulating reactant water removal technique. 


MODEL DEVELOPMENT 


A schematic diagram of the cell for which the model was developed is shown in fig- 
ure 1. The cell shown in this diagram is a matrix cell of the same type used in an ex 
perimental water- rejection dynamics test program previously conducted (ref. 2). In 
this cell the electrodes consist of a nickel screen coated with a mixture of platinum 
black catalyst and Teflon binder, and the electrolyte is an aqueous potassium hydroxide 
solution. The product water is evaporated from the electrolyte contained in the hydro - 
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Figure 1. - Fuel cell. 


gen electrode and removed from the cell by flowing a controlled excess of humid hydro- 
gen through the hydrogen supply chamber. 

The mass transport model that was developed is one dimensional. This character- 
istic is implicit in the macroscopic or lumped parameter analysis that was used in the 
derivation. 

The model is analytical in all aspects except in the simulation of the effect that cell 
operating temperature has on the water- rejection process. In lieu of deriving a com- 
plete heat- transfer model to provide a totally analytical simulation, the thermal charac- 
teristics and their related effects on the water- rejection process were synthesized em- 
pirically when testing the validity of the analytical model on the analog computer; that 
is, for any particular perturbation, experimental cell- temperature response data were 
used to obtain a corresponding correction in electrolyte vapor pressure. Since the rate 
of water rejection is a function of the electrolyte vapor pressure, the thermal effects 
were accounted for by superimposing the temperature correction on the program- 
calculated vapor pressure. 

Accounting for the temperature effect in this manner assumes that temperature 
gradients within the cell are negligible. This assumption was justified by experimental 
thermal response data taken during various load- step transients in which a negligible 
difference was observed between the temperatures of the reactant- side surfaces of the 
two electrodes. 

The zones selected for deriving the model are depicted in schematic form in fig- 
ure 2. The transport of water and OH ions that occurs within the electrode- matrix 
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structure is included in the analysis, as is the rejection of water vapor into the hydrogen 
supply stream. The transport of water within the reactant chamber is considered to be 
at steady- state because the lag or retention time of water vapor within the chamber is 
negligible. 

To represent the electrodes, a flooded agglomerate configuration is used (ref. 3). 
The theory proposed in the development of this electrode model states that the catalyst 
present in Teflon bonded platinum electrodes tends to form porous agglomerates that 
flood with electrolyte. The presence of Teflon creates hydrophobic gas filled regions 
that surround the agglomerates and provide passages for reactants to diffuse into the 
available catalyst reaction sites. In the case of the hydrogen electrode, these regions 
also provide a path for the evaporated product water to diffuse into the flowing hydrogen 
stream. To approximate the electrode geometry, the agglomerates can be considered 
to be electrolyte- filled cylinders that extend from the matrix to the gas- side surface of 
the electrode (see fig. 2). The electrolyte completely fills the matrix plus a portion of 
each electrode, thereby flooding not only the agglomerates, but also a portion of the void 
volume surrounding the agglomerates (subsequently referred to as the gas- filled hydro- 
phobic regions, GFHR). Changes in electrolyte volume due to changes in electrolyte 
concentration are taken up by the gas- filled hydrophobic regions; the result being a var- 
iation in the flooded volume of each electrode. It is assumed that under all conditions 
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both electrodes are flooded by the same amount; that is, the ratio of filled to unfilled 
volume is equal. 

Within the cell configuration depicted in figure 2, the four mass transport zones are 

(1) The hydrogen supply chamber: this is the rectangularly shaped gas passage 
through which the humidified hydrogen stream is supplied to the cell 

(2) The bulk electrolyte zone: this is composed of the volume of KOH- filled matrix 
plus the volume occupied by the KOH that floods into the originally gas- filled regions 
that surround the flooded electrode- agglomerates 

(3) and (4) The hydrogen and oxygen reaction zones: these consist of the Siam total 
of the volume occupied by the agglomerate cylinders that extend from the electrolyte- 
filled portion to the gas- side surface of each electrode. 

The mathematical analysis is based on the continuity of water and hydroxide ions 
present in the zones described. The resulting mathematical simulation, the derivation 
of which is presented in detail in appendix A, is expressed by five nonlinear first-order 
differential equations. 

The model equations, in the form in which they were programmed for solution on an 
analog computer, are the following continuity equations: 

(1) Water in the hydrogen reaction zone 


' m h 2 o-h 2 m h 2 o-b n 


dM 


H 2°- H 2 

dt 


- 1.87x10 (I) - CA (P* - P g }i ) - d h 9 0-K0H A Z 


Hr 


V 


B 


(2) Water in the bulk electrolyte zone 


(A20) 


/ M H n O 


dM 


h 2 o-b _ 


dt ■ d h 2 o-koh a z 


2~- H 2 Mh 2°- B 


'H 0 


J B 


- D 


HgO-KOH^Z 


V 


m h 2 o-b m h 2 o-o 2 


A 


(A22) 
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(3) Water in the oxygen reaction zone 


dM H 2 0-0 2 

dt 


- d h 2 o-koh a z 


( U h 2 o-b m h 2 0-o 2 \ 



0. 935X10“ 7 (I) 


(4) The hydroxide ions (OH“ ) in the hydrogen reaction zone 


(A24) 


dM OH"-H 2 

dt 


/“OH--B M OH--H 2 \ 



1. 78X10" 7 (I) 


(5) The hydroxide ions (OH”) in the bulk electrolyte zone 


(A28) 



(A29) 


The symbols used in the preceding equations are defined in appendix B, and the values 
of the various cell parameters are presented in tabular form in the model derivation of 
appendix A. 

The dependent variables in the model equations are the masses of water and OH” 
present in the zones for which the respective continuity equations were written, and the 
independent variable is time. The mass flow expressions in the continuity equations are 
provided by applying Pick's law to the various model zones since the transport within the 
cell is diffusive in nature. The driving functions are the load current I and the partial 
pressure of the water vapor in the inlet- hydrogen stream P . . The load current is a 
measure of the rate of water production, and the inlet stream vapor pressure is a meas- 
ure of the water content or degree of humidification of the hydrogen stream. 

Nonlinearities are present in the model because various cell parameters, which 
appear in the coefficients of these equations, are functions of the system dependent var- 
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iables, which, in this case, axe the transport zone constituent masses; that is, param- 
eters such as effective diffusivity, diffusion path length, and mass-transfer area are 
functions of the instantaneous electrolyte volume, which, in turn, is defined by the in- 
stantaneous mass makeup of the electrolyte in the various transport zones. 


MODEL VERIFICATION 


Verification of the mathematical model was accomplished by analytically reproduc- 
ing experimentally obtained water- rejection response data. The verification consisted 
of comparing experimental data with the change in the rate of water rejection as meas- 
ured by the term 


( C ) A Z-H 2 ( P * 




which appears in equation (A20). This expression actually represents the rate at which 
water is rejected from the hydrogen electrode into the hydrogen stream. However, by 
applying the assumption that no transport lag occurs in the gas chamber, this term can 
be used as a measure of the change in the flow rate of water exiting the cell and can then 
be compared with the water- rejection data obtained in the experimental study (ref. 2). 

Both experimental and analytical data for a typical water- rejection response to a 
step in the rate of water production (load current) are shown in figure 3. In this figure 
the model correlates the experimental data to within 5 percent. It predicts the correct 
peaking and settling time of what is basically an underdamped second- order response, 
and it also approximates the slight ripple or higher- than- second- order characteristics 
that is superimposed on the basic second- order curve. 

A typical model and experimental response correlation for a step in inlet stream 
humidity is shown in figure 4. Again, as for the power- step case, the maximum devia- 
tion between experimental and analytical data is less than 5 percent. 

The changes in the experimentally observed responses caused by varying the ampli- 
tude of the input perturbation as well as the changes brought about by variations in pre- 
transient operating conditions are also correctly predicted by the model. 

The effect that the amplitude of the input step has on the water- rejection response 
is shown in figure 5. The curves of this figure are plots of normalized peak outlet water 
flow against load- step magnitude. Normalization was accomplished by dividing peak 
outlet water flow values by final (steady state) values. This was done for load current 
steps of several different magnitudes. The experimental curve, which was taken from 
the report on the matrix cell work (ref. 2), is compared with the corresponding analyti- 
cal version obtained with the model. 
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Figure 3. - Comparison of model and experimental outlet-stream water flow response to step in load 
current. Operating conditions: current step, 20 to 60 amperes; inlet stream hydrogen content 
0.05 kilogram per hour (0. 12 lb/ hr); water content of inlet stream, 0.05 kilogram per hour 
(0. 12 Ib/hr). 
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Figure 4. - Comparison of model and experimental outlet-stream water flow response to step in water content of 
the inlet hydrogen stream. Operating conditions: inlet-stream water content step, 0. 082 to 0. 11 kilogram per 
hour (0. 18 to 0. 24 lb/hr) ; inlet stream hydrogen content 0.05 kilogram per hour (0. 12 lb/hr) ; load current. 

41 amperes. 



Figure 5. - Normalized peak outlet water flow rate as function of mag- 
nitude of load current step transient. Operating conditions: inlet- 
stream water content 0.082 kilogram per hour (0. 18 lb/hr); inlet- 
stream hydrogen content, 0.05 kilogram per hour (0. 12 lb/hr); 
initial load current, 20 amperes. 
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The plot of this peak- to- final ratio against current step magnitude identifies the 
water-rejection process for steps in cell current as being nonlinear (a linear process 
would result in a straight line). The close agreement between the curves of figure 5 
indicates that the model correctly predicts this system nonlinearity. 

Typical of the changes in response caused by variations in initial or prestep operat- 
ing conditions are the data shown in figure 6(a). The two curves in this figure are re- 
sponses to the same size power step taken at two different prestep inlet stream humid- 
ities. Since the prestep values and final steady- state water outflow values for the two 
differ from each other, the data are plotted on a normalized or percent- change basis to 
allow a direct comparison and to illustrate the response difference between the curves. 
The curves in figure 6(a) show that, as the inlet humidity is increased, the response 
overshoot decreases. The model simulates this same effect as seen from the data 
plotted in figure 6(b). 

This response deviation is revealed by the model to be the result of system nonlin- 
earities that stem from two matrix cell characteristics. These characteristics, which 
are discussed at greater length in appendix A, are the inlet hydrogen stream dependency 
of the electrode agglomerate mass-transfer area and the relative amount of electrolyte 
in the electrodes and matrix (see fig. 2). In the model these characteristics are trans- 
lated into coefficients that are functions of the dependent variables and produce the non- 
linear results. 

To determine whether this water- transport simulation is general enough to predict 
the basic response characteristics of a Bacon cell, the model was adjusted to describe 
the Bacon cell used in the earlier water- rejection study (ref. 1), and several of the ex- 
perimental perturbations were simulated. The adjustments that were made entailed 
changes in operating conditions, geometric parameters, and defining cell characteristics. 
The transport processes that are modeled were not changed, and thus the basic simula- 
tion remained intact. 

An example of the type of cell characteristics taken into account is the difference 
between the "simple pore" configuration that represents the dual- porosity, sintered- 
nickel Bacon cell electrodes and the "flooded agglomerate" concept that simulates 
the Teflon bonded platinum catalyst electrodes of the matrix cell. In the model this 
difference was accounted for in the hydrogen reaction zone mass- transfer area. The 
mass-transfer area for the Bacon cell is merely the surface area of the electrode, while 
for the matrix cell it is the sum total of the cylindrical surface areas of the agglomer- 
ates in the nonflooded portion of the electrode. 

A typical Bacon cell experimental response to a step in load current or the rate of 
water production along with the corresponding model response curve are shown in fig- 
ure 7. Again, as for the matrix cell, the experimental and analytical response data are 
in good agreement. 
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Figure 7. - Comparison of model and experimental Bacon cell outlet stream water flow response to a step in load 
current. Operating conditions: current step, 10 to 15 amperes; inlet stream dry hydrogen at 0.0058 kilogram 
per hour (0.0128 Ib/hr). 


Thus, since the model reproduces the basic response curves of both types of exper- 
imental transients for both cells plus all the various operating effects and system non- 
linearities, it is concluded that it is a good generalized approximation of the water trans- 
port in Bacon and matrix cells which employ circulating reactant- stream water removal. 


DISCUSSION 


The basic response characteristics for both types of experimental steps for both 
cells were reproduced analytically as were all the various operating effects and system 
nonlinearities. Since this was accomplished, the validity of the model is established 
even though it was necessary to make approximations in sizing various undefined cell 
parameters such as the electrode and matrix tortuosity, the OH" ion transport coefficient, 
and the radii of the cylindrical electrode agglomerates. The "curve fitting" latitude 
provided by these approximations might be sufficient to reproduce a single transient; 
however, it is obviously not sufficient to reproduce all the response characteristics. 

An interesting matrix cell characteristic predicted by the model is that the water- 
rejection overshoot or drying evidenced by the load- step data cannot be totally attributed 
to a thermal effect; that is, the effect in which a cell temperature rise produced by a 
load step increases the electrolyte vapor pressure, which, in turn, causes water to be 
driven from the cell. Rather, the drying that occurs is, in part, the result of an inher- 
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Figure 8. - Model outlet-stream water flow response to step in load current with cell operating isother- 
mally. Operating conditions: current step, 20 to 60 amperes; inlet stream hydrogen content 0. 05 
kilogram per hour <0. 12 lb/ hr); water content of inlet stream, kilogram per hour (0. 12 Ib/hr). 


ent characteristic of the cell. This is illustrated in figure 8 which is the model response 
for the load step shown in figure 3 with the temperature effects left out. The overshoot 
or drying effect, though not of the correct magnitude, is still present. This indicates 
that the basic transport process in a matrix cell is a higher- than- first- order phenome- 
non, an important fact when determining the characteristics of a system to control the 
water inventory for this type of cell. 

A possible physical explanation is that the higher- than- first- order characteristic 
results from an imbalance in the diffusive resistance that exists across the cell struc- 
ture. A step increase in load current produces an increase in the rate of water genera- 
tion at the anode and an increase, one- half as large, in the rate of consumption at the 
cathode. If, as in the model, the resistance to water transport through the bulk electro- 
lyte toward the cathode is higher than it is from the hydrogen reaction zone into the gas 
stream, the cell would have a tendency to reject more than one- half of the amount being 
generated at the anode. To satisfy the increased cathodic demand, water would be with- 
drawn from the bulk reservoir; the end result being on the reestablishment of equili- 
brium, a net removal of water from the cell. 

In regards to generality or applicability to other types of cells, it is believed that 
the model could be used as a basis from which various additions or deletions would be 
made to account for differences between the type of cells modeled and the particular cell 
in question. Such additions to the existing model could represent such characteristics as 
a sizeable thermal gradient caused by liquid cooling, or the concentration gradients that 
accompany a static water removal system. 
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CONCLUSIONS 


From the development of the mathematical water- transport model the following 
conclusions can be drawn: 

1. A generalized mathematical simulation was developed for the water transport in 
Bacon and matrix cells that employ circulating reactant stream water removal. The 
validity of the model was established by the fact that it analytically reproduces both the 
linear and nonlinear transient response characteristics of both types of cells. 

2. A matrix cell characteristic predicted by the model is that the basic transport, 
exclusive of operating temperature effects, is a higher-than- first- order process. The 
significance of this is that the electrolyte dr ying that accompanies a matrix cell load 
increase can be partially attributed to this higher-order effect and would thus be present 
to some degree even if the cell temperature were held constant during transient opera- 
tion. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, September 27, 1971, 

113-34. 
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APPENDIX A 


DERIVATION OF THE MATHEMATICAL MODEL 

The electrochemical reaction at the anode 

H 2 + 20H” - 2H 2 0 + 2e" 

dictates that for every mole of hydrogen reactant consumed by the cell, two moles of 
water are formed. At steady state, one mole of product water must be rejected from 
the cell, and one mole is transported across the cell to be consumed in the overall 
cathodic reaction which is 

°2 

— + H 9 0 + 2e - 20H 
2 z 

The hydroxide ions, which are the product of this reaction, are, in turn, trans- 
ported to the hydrogen electrode for consumption in the anode reaction. 

To simulate the cell mass transport that occurs during steady- state and transient 
operation, the continuity (mass conservation) equation was applied, on a lumped param- 
eter basis, to the hydrogen, oxygen, and bulk transport zones defined in figure 2. 


CONTINUITY EQUATION FOR WATER IN THE HYDROGEN REACTION ZONE 


The continuity equation for water in the hydrogen reaction zone volume states that 
the rate of mass accumulation equals the rate of mass generation minus the rate of out- 
flow. The equation can be expressed as 


dM 


h 2 o-h 2 

dt 


= w- 


R 


w„ 


w 


E-H 0 


(Al) 


where w^ represents the generation rate w g the mass outflow from the reaction zone 

into the hydrogen stream, and w^jj the outflow into the bulk electrolyte zone. 

2 

Inherent in a zonal or lumped parameter analysis such as this is the assumption 
that there are no gradients in the variables or cell parameters of any transport zone. 
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Rather, each parameter or variable within a zone is defined by an average value which 
is assumed to extend throughout that zone. In terms of the analysis, this states that the 
transport dynamic effects internal to the various zones are neglected. 


DERIVATION OF RATE EXPRESSION FOR FLOW OF 
WATER VAPOR INTO HYDROGEN STREAM 


The first of the rate terms to be derived will be w„, since it depends to a higher 

S 

degree on the details of the particular cell under consideration than do and Wg_ jj 

2 

This sensitivity to cell detail plus the fact that it is necessary to describe w g quite 
accurately to obtain satisfactory agreement with experiment make it impossible to de- 
rive a single, completely general expression for w g . It will be shown, however, that 
cells of quite different character can, in many cases, be described by the same expres- 
sion for w_. In any case, a minimum amount of response data is needed to decide 
s 

upon a specific form of w g . 

In the next section two expressions for Njj q will be considered in the light of the 

2 

flooded agglomerate model geometry, where q is the (local) flow rate of K 2 0 vapor 

2 

out of the electrode per unit electrode geometric area. The first expression pertains if 
transport of water vapor in the gas-filled hydrophobic regions (GFHR) of the anode is 
controlled by a resistance associated with the liquid-gas interface. The second expres- 
sion pertains if vapor phase diffusion in the z direction (perpendicular to the electrode 
surface) controls transport. Both expressions turn out to have the form 


N, 


h 2 ° 


= h(P* - P g ) 


where P* is the local (x- dependent) electrolyte equilibrium vapor pressure and P_ is 
the corresponding local HgO partial pressure in the bulk gas phase in the gas cavity. 

Note that P* is the mean vapor pressure over the z- direction. The parameters 
appearing in h, the local mass-transfer coefficient, differ according to which resistance 
is controlling. 

* 

Next, three possible distributions of N^ q and P over the face of the electrode 

will be considered. The first represents one extreme - that of constant (i. e. , position 

independent) Njj q, where P* is constant and the gas in the gas cavity is perfectly 

2 ^ 
mixed. The second case is a condition of maximum variation of q. Here, too, P 

2 

is constant, but the bulk gas is in plug flow. The second case may be contrasted with 
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the third, wherein the bulk gas is in plug flow but the vapor pressure P* varies with 
the stream partial pressure P g such that d/dx (P* - P g ) = 0, where x is the flow- 
path directed distance from the cell gas inlet. 

All three cases yield expressions of the form 

W -i 

•sr- = xf/ (u) — (P* - P .) 

W T n S>1 

where x}/ (u) is a dimensionless overall mass- transfer coefficient, n is the mean total 
pressure in the gas cavity, P* is the overall mean equilibrium vapor pressure, P c . 
is the partial pressure of HgO in the gas stream entering the cell, and w T is the total 
molar flow of gas into the cell (note that strictly speaking the forms of w„ to be derived 

o 

apply only to cells in which product water is removed by excess hydrogen flow. So- 
called static water removal methods may, of course, be treated similarly but would 
yield somewhat different forms for w 0 ). 

The dimensionless quantity u is given by 


hA z 



n 


where A z is the geometric electrode area. The three different cases mentioned pre- 
viously yield different functional dependencies of \f/ on u. 

Finally, in the section Form of the Rate Expression the similarities and differences 
in the various resulting expressions will be discussed briefly, as will be the process by 
which one may arrive at a form applicable in a given situation. The form adopted for 
the simulation of the data given in reference 2 will be used as an example. 


Water Vapor Transport In the Gas-Filled Hydrophobic Regions 

Conservation of HgO vapor in a differential volume of the GFHR yields 


‘“V N 9 

dZ’ 1 cp 


(A2) 


where the accumulation of HgO vapor and net transport in directions other than perpen- 
dicular to the electrode surface (the Z’- direction) is neglected. In equation (A2) 
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Nr o * s m °l ar fr ux °f HgO in the gas -filled hydrophobic region in the Z’ direction - 
2 

the direction parallel to the axes of the agglomerates and towards the gas cavity; ISh is 
the molar rate of HgO transport across a unit area of gas- liquid interface or agglom- 
erate surface; 6 is the interfacial area per unit electrode geometric volume; and (p is 
the gas- filled cross sectional area in the Z ’-direction per unit electrode geometric sur- 
face area. Both 6 and (p may be assumed to be independent of Z in accord with the 
flooded agglomerate model configuration. 

Interfacial resistance controlling . - The chief transport resistance is postulated to 
be associated with the gas liquid interface. From this point of view the following 
assumptions are made: 

(1) There are no Z’-direction HgO concentration gradients in the liquid or gas phase. 

Thus, = P g for all Z', where P g is the partial pressure of HgO in the gas 

cavity immediately above the portion of electrode under consideration. 

(2) The molar flux N. is given by 


N i =k i 


P* - P H2 o(Z’) 


= kj(p* - P s ) 


where k- is the mass- transfer coefficient associated with the local gas- liquid interface 
. and P* is the equilibrium HgO vapor pressure associated with the local liquid phase. 

(3) All HgO evaporates from the lateral surfaces of the agglomerates. Thus, 

NJj q = 0 at Z’ = 0 (the base of the agglomerates). 

Introducing these assumptions into equation (A2) yields 

dN H 2 0 e 

=— k.(P* - PJ 

dZ ' cp 1 s 

Integrating and applying assumption (3) yield 

N H 2 0-f k l< P *- p s> Z ' 

The rate of HgO flow into the gas cavity per unit electrode geometric area q is 
given by 2 


n h 2 o _ ^ n h 2 o 


= 0k i Z H ( P * - P s) 


(A3) 


Z’=Z 


H 
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where Zg is the nonflooded agglomerate height. 

Vapor diffusion controlling . - Here the transport resistance is postulated to be the 
Z-direction diff usion of water vapor. The following assumptions are made: 

(1) The molar flux is constant over the entire gas- liquid interface area; that is, 
N- is independent of Z*. 

(2) The molar flux N' „ results from equimolar countercurrent diffusion (i. e. , 

h 2 u 

there is an equal but oppositely directed flux of hydrogen). Thus, 


n h 2 o 


Dg P H 2 0< Z > 
RT dZ’ 


(A4) 


where is the binary diffusion coefficient for HgO - H 2 , R is the gas constant, and 

T the absolute temperature. 

(3) The molar flux N-^ q is zero at Z' = 0. 

2 


Applying these assumptions to equation (A2) and integrating twice yield (recalling 


that P R Q (Z) 
2 


P g at Z’ = Z„) 


P H*0< Z > “ P s 


N.0RT / 9 

— — (zi 

2 <PD V 



(A5) 


Since in this approach the controlling transport resistance is Z-direction diffusion, 
the HgO vapor in the GFHR may be considered to be in equilibrium with the surrounding 
liquid. Thus the mean partial pressure of HgO vapor in the GFHR may be identified with 
the mean vapor pressure of the liquid in the local agglomerates. Mathematically, 

/* Z H 

P * = P H,0 < Z) ‘7- / P H,0< Z > dZ ' 

2 Z H J 0 1 


Thus, 


P * = P s + 


N.0RT 


3 <pD 


J H 


g 


This result may be solved for N^, which then may be eliminated from equation (A5). 

The flux Njj q may now be obtained from diffentiating the resulting expression: 
2 
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n h 2 o - ^ n h 2 o 


Z’=Z 


H 


<pD g dP H 2 0 
RT dZ’ 


Z’=Z 


H 


Thus 


N 


H 2° RTZ 


3D (p 

— (P* - P s ) 


(A6) 


H 


Alternate expressions are now available for q, the assumption of interface con- 
trolled transport leading to equation (A3) and the assumption of diffusion limited trans- 
port yielding equation (A6). Note that the driving force is the same but that the form of 
the coefficient differs. To summarize 


n h 2 o - h < P * 


P s> 


(A7) 


with 


h = h., = 


j = 0kjZjj for interface control 


or 


hHh 2 = 


3D 


RTZ 


H 


for diffusion control 


where in either case h may be thought of as a local mass-transfer coefficient. 


Distributions of Local Mass Flux and Vapor Pressure 


In the preceding sections, care was taken to use the modifier "local" when dis- 

jjf 

cussing the quantities appearing in the expressions for Q . Indeed, q and P 

may vary from place to place on the electrode, depending in a complicated way on such 

things as current density, temperature, and gas cavity partial pressure distributions. 

It is beyond the scope of this study to consider such dependency in detail. Rather, three 

representative cases will be postulated, which represent the spectrum of possible varia- 

* * 

tion in Njj q and P : (a) Njj q and P independent of position, (b) maximum varia- 
2 2 
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tion in q, and (c) maximum variation in P 
2 


Case (a) may be associated with a condition of maximum or '’perfect” mixing in the 
bulk gas phase in the gas cavity. "Perfect mixing" is a commonly used term referring 
to a condition where there are no concentration gradients and where the outlet stream 
or streams are identical in composition to the gas in the gas cavity. 

Cases (b) and (c), on the other hand, may be associated with a condition of minimum 
mixing or plug flow in the gas cavity. The term "plug flow” refers to a situation in 
which there is a definite direction (e.g. , the X direction) to the flow and in which 
there is no mixing in the flow direction. Plug flow, then, yields the largest possible 
partial pressure gradients in the gas cavity. (For further discussion and examples of 
use of these mixing models see any text on chemical reactor design (e.g. , ref. 4). ) 

The type of flow actually present in a given cell depends, of course, on the design 
of the gas cavity. Thus in some cases, it will be clear which type of flow exists. But 
even when the type of flow is not clear, one of the cases considered should be sufficient 
to describe the dynamic behavior of the cell since it will be seen that the resulting ex- 


pressions for w c 


And under some conditions 


w_ will be the same 

fc> 


, g are quite similar 
irrespective of the type of flow. 

Case a: flux and vapor pressure i n dependent of posi t ion (perfect mixing) . 

and N tt is 


The 


general relation between w, 


s 


■h 2 ° 






n h 2 o dA 


(A8) 


where dA is an increment of electrode surface area and A^. is the entire electrode 
surface area. Assuming that h in equation (A7) is independent of position and incor- 
porating the perfect mixing characteristic (constant P„), equation (A8) may be inte- 
grated in the perfect mixing case to give 

w g - hA z (P* - P g ) (A9) 

where 

P* = JL f P* dA 
A z A z 

that is, P* is the overall mean vapor pressure associated with the reaction zone - the 
quantity needed for the complete lumped parameter model. 
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It remains to eliminate P„ from the expression for w„. This can be accomplished 

b b 

by means of the equation for conservation of H 2 0 vapor in the gas cavity- Neglecting 
the accumulation term results in 




W„ = W 


:out 

h 2 o 


where the 


w h 2 0 


’s are molar flow rates. 


(It may be shown, with an order- of- magnitude 


analysis, that the accumulation term is generally negligible because of the very small 
residence times associated with the flowing stream. ) 

The molar flows may be related to partial pressures by the relations 




n 





and 


w 


out 

HoO 


s ~ 

— Wn 


n 


where w,p is the total molar flow rate through the cell and II is the total pressure. 
Both w T and n generally vary only slightly from inlet to outlet, so w T /n may be 
assumed constant. Note that the expression for Wjj q includes the perfect mixing 

characteristic that P^ U q = P g . 

2 

Combining these relations with the expression for w„ yields 

s 


^w T + hA 7 (P* - P_) =— w n 


n 


n 


This may be solved for P : 

b 


Wr 


_ p* n_ 


W r 


n 


+ hAr 
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And this expression may be substituted into equation (A9) to yield 





n 


u 

1 + u 


(P* - p 



(A10) 


where 


u = 


hA z 



n 


Note that equation (A10) is of the form 


w_ 


s 


w 


_T 

n 


- i// (u)(P* - P s> .) 


where ^ may be thought of as a dimensionless overall mass-transfer coefficient. 

Case b: maximum variation in flux (plug flow with constant vapor pressure) . - Plug 
flow is characterized by the fact that no element of the fluid overtakes another element; 
that is, there is no mixing of the fluid in the direction of flow, although there may be 
mixing in the lateral direction. 

If flow is in the X direction, a HgO vapor mass balance on a differential volume 
element of the gas cavity yields 


dw 


h 2 ° 


dX 


= n h 2 o* 


where Wjj q is the molar flow of HgO vapor and l is flow path width. Note that the 
2 

accumulation term is neglected. 

The molar flow may be related to the partial pressure by 


~ w T 

Wtt rw = P c 

h 2° s n 


where w T /n may again be assumed constant. 
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Substituting this relation and the general form of Nr Q (eq. (A 7)) into the mass 
balance yields 


dP 

s 

dX 


4P (P* - P s ) 

W™ 


(All) 


This relation is characteristic of the plug flow condition. 

* 

The position dependency of P must now be introduced into equation (All). It may be 

shown that the variation in Nr q is maximized, subject to the constraints that Nr q 

2 2 
and dP* /dX be zero or positive, if P* is constant (these constraints arise from 

practical fuel cell operating requirements). Thus, for maximum variation in Nr q 

with position, the solution to equation (All) is ^ 



which employs the boundary condition that P„ = P 0 . at X = 0. 

S Sj 1 

Substitution of this expression for P g into the expression for Nr q (eq. (A7)) gives 
Nr o with maximum position dependency. ^ 

The general relation between w g and Nr q (eq. (A8)) may now be integrated for 

2 

this particular form of Nr q with the result that 

2 


where, again, 


(V 

\v 


s 


= (1- e" u )(P" - P s>i ) 


w 


_T 

n 


hA 



n 


and the expression for w„ is of the form 

o 


(A12) 
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^r-=*( u)(P* - P s i) 

w T S)1 

n 


Case c: maximum variation in vapor pressure (plug flow with constant flux) . - In 

this case the flux N H q is to remain constant, and the flow in the gas cavity is plug 

2 ^ 

flow. Equivalently, d(P - P )/dx is zero. 

s 

Equation (All) again pertains; in this case the right hand side is constant, and so 
integration yields 


P s =P s i + ^ < p * - p s ) x 
s Sjl w T s 


At the same time equation (A8) may be integrated for constant q to give 


S s = n h 2 0 A Z = ' Lh <P* - P s ) (A13) 

where L is the flow path length. 

Combining these two expressions yields 


nw x 

P = P . + _ ii 

S S;1 Wrr. L 


Since it is desired to express w g in terms of P , the constancy of Njj q in this 
case may be taken advantage of to write equation (A13) as ''' 


w_ = Ar 7 h(P* - P Q ) 


The expression for P„ may be integrated to yield P : 

s s 


- 1 nw s 

p„ = p„ . + - 


s S ,1 2 


Wn 


These two expressions may be combined to eliminate P„ and may be solved for 

s 


w g to give 
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s _ u /^* 


W T 1+S 


(P* - P s?i ) 


(A14) 


where once again 


and the 1 relation again is of the form 


^=*(u)(P* - P s i) 


Form of the Rate Expression 

It has been seen that all the cases considered in the preceding section yield the same 
type of expression for w g , namely, 

w __ 

i//(u)(P* - P .) 

1 


From the section Water Vapor Transport in Gas-Filled Hydrophobic Region, the 
alternatives for h are 


h s hj = 0k^Zjj for interface control 
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for diffusion control 


h s 


* RTZ H 


And from the section Distributions of Local Mass Flux and Vapor Pressure, the three 
alternatives for ip (u) are 


1p(u) S 1p t =—±— 

(A10) 

1 + u 


ip(n) = ip 2 = 1 - e" u 

(A12) 

ip ( 11 ) = ip 3 =— — 
1 . u 

(A14) 


2 


Figure 9 shows the behavior of ip^, ip 2 , anc ^ ^3 over several decades of u 
values. 

The first property common to the ip's, regardless of functional form, is that they 
are monotonic with u and, hence, with h. This makes it possible to determine from 
a small amount of experimental data which type of transport takes place in the GFHR, 
that is, which form of h (hj or hg) is applicable to a particular fuel cell. 
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This determination may be made as follows: Two experimental transients are re- 
quired, of the sort shown in figure 6. (Based on data from ref. 2. ) These data in effect 
illustrate differences in w g /w T due to different initial conditions - the high inlet stream 
humidity corresponding to a relatively dilute electrolyte or ''wet" condition and the low 
humidity to a concentrated electrolyte or "dry" condition. The wet and dry conditions, 
in turn, may be identified with smaller and larger reaction zone thicknesses, or agglom- 
erate heights (Zjj), respectively. Since h^ is proportional to Z H while 1^ is pro- 
portional to l/Zjp if hj is the correct form for the cell in question, then w g should 
be generally greater for the transient with the dryer initial condition. Conversely, if 
hg behavior is controlling, w g should be greater during the transient having the wetter 
initial condition. The sense of the word "greater" is evident from figure 6, which indi- 
cates that the cell of reference 2 behaves in the fashion described by h^; that is, the 
controlling GFHR transport resistance is associated with the gas- liquid interface. The 
hj form for h was, in fact, used in all simulations of reference 2 data. 

The quantities in hg may a priori be assigned numerical values (the calculation of 
geometric parameters associated with the flooded agglomerates will be discussed later), 
but the interfacial mass-transfer coefficient that appears in h^ is difficult to esti- 
mate from theory (although it should not vary during a transient). Because of this fact 
it was necessary, when simulating the data of reference 2, to use the analog computer 
simulation to calculate the initial steady- state value of xp (u). This was accomplished by 
putting the computer into operation without introducing a perturbation function into the 
program. Under these conditions the computer was made to calculate the pretransient 
initial steady- state value of xp, which, in turn, fixed the value of k. since k. is con- 
stant throughout a transient. 

A second feature of figure 9 is that xp xp^, and xp^ become identical for small 
values of u. More particularly for u < 10” xp^ = xp 2 = ^3 — u. Thus, for a cell with 
poor GFHR transport and a high gas- stream flow rate w_ becomes independent of de- 

gree of mixing in the gas cavity and of vapor pressure gradients. 

9 

It is also obvious that for u > 10 ^(u) becomes independent of u. The dynamic 
behavior of the cell thus becomes linear if good GFHR transport is coupled with low gas 
flow rates. It can be seen that this corresponds to equilibrium between the exiting gas 
and the liquid at the exit. 

A final feature of figure 9 is that the x p dependency on u is readily linearized over 
small ranges of u. It is easily shown that over any range of u where 

( u max - u min^ * 

the following linear approximations may be made: 
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I 




^ + u initial. 


^9 ~ e 


initial u 


u initial 


This linearization simplifies calculations considerably, since the technique may be 
applied - that of letting the programmed model calculate xp values based on initial con- 
ditions. It was, in fact, found that the linearized form could be used to simulate all the 
data from reference 2. Note that, if xp may be linearized and if the model is pro- 
grammed to calculate initial condition ^'s, then it is unnecessary to decide which form 
of xp pertains; that is, details of gas cavity flow and vapor pressure distribution be- 
come unnecessary. 

The following expression, then, may hold in many cases, regardless of details on 
cell design: 


w g = c 


h- o?*-p M > 


where c is a coefficient dependent on initial conditions only and all nonlinear behavior 
of w g is centered in h, which, in turn, reflects changing agglomerate geometry. 
TaMng advantage of the relation 


A Z Z H 


allows u to be written for h s h^ as 


A Z-H 2 k i 
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Here, again, C is a coefficient generated by the programmed model for each pretran- 
sient initial condition. 


Calculation of Geometric Parameters 

The agglomerate surface area is given by 

A Z-H 2 - A Z Nnr o + A Z N2nr o ZH 

2 

in which A z NIIr 0 represents the total base area of the agglomerate cylinders abutting 
the gas chamber and A z N2nr Q Zjj the total lateral area. The number of agglomerates 
per square centimeter of electrode N was calculated from the agglomerate model de- 
finition of macroporosity (percentage of electrode not occupied by agglomerates). The 
equations used to define this void percentage are 

<p= 1 2 

pH(l - e) 

and 

(p = 1 - NUr 2 

7 

By equating these two expressions and solving for N, a value of 1. 64X10 cylinders per 
square centimeter was determined, and with this the area was calculated to be 

A z _ h = (1.7 + 1300 Z H )cm 2 
2 

In the programmed model, the calculation of the agglomerate height, Z H is based 
on a knowledge of the water and potassium hydroxide masses in the bulk transport zone 
at any instant of time. The model equations provide continuous calculations of the water 
mass and the OH” component of the bulk zone electrolyte. To determine from the OH" 
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mass calculation the amount of KOH present at any time, the assumption of electrolyte 
electroneutrality (ref. 5) was employed. This approximation is based on the fact that 
the ionic species contained in an electrolytic solution, such as that in the cell being 
modeled, have sufficient mobility to effectively neutralize the charge at any point within 
the solution. In terms of the model this defines the molar concentration of K + and OH" 
of any transport zone to be equal. Therefore, for the bulk transport zone 


M 

m oh--b k + -b 

17V b 39V b 


(A16) 


and thus the bulk zone KOH mass is 

M KOH-B = M k +_ b + M OH“-B = + X ) M OlT -B (A17) 

The water and KOH mass values are then used to determine the bulk zone weight 
fraction Xj^qh an d hence the bulk mixture density from the empirical relation 

P E = 0. 950 + 1. 045 X K q H 

The data used in deriving this expression are presented in reference 6. The bulk zone 
density and liquid- mass composition define the bulk liquid volume from which the bulk 
zone thickness is calculated to be 


Z B = 


A Z P 


where P is the bulk zone mean porosity. By subtracting Zg from the total thickness 
of the electrode-matrix structure, the sum of the heights in the nonflooded portions of 
the two electrodes is obtained. The nonflooded portion of each electrode Zg is one- 
half of this sum, since the bulk electrolyte volume that extends into the electrodes is 
assumed to be distributed equally between them. 

Thus, through the agglomerate- height, electrolyte-volume relation, the area calcu- 
lation becomes functionally related to the component water and KOH masses of the elec- 
trolyte. In terms of the analysis this represents an equation coefficient A ZB , which 

is a function of the dependent variables of the describing equations; hence a nonlinearity 
is introduced. In the model this accounts for part of the nonlinear characteristic that 
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was produced experimentally by varying the preload- step- transient inlet stream humidity 
(fig. o); that is, a variation in inlet stream vapor pressure (inlet hydrogen stream 
humidity) produces a change in the mass composition or electrolyte volume, which, in 
turn, is reflected into the value of the area Ag_ jj • 


DERIVATION OF RATE EXPRESSION FOR DIFFUSION 
OF WATER INTO THE BULK LIQUID 


In the general hydrogen reaction zone continuity equation (ea. (Al)), the expression 

defining the flow of water from the reaction zone into the bulk zone Wg_ ^ is defined 

2 


by Fick’s law. However, Fick’s law in this case is applied to mass transport in the 
liquid state. The equation is 


dp H 2 0 


-»E-H 2 = ^O-KOI^Z — 


On a lumped parameter basis this is expressed as 


/ m h 9 o 


w e-h 2 _ d h 2 o-koh a z 


2 w-H 2 m h 2 o-b 


Hr 


J B 


) 


(A18) 


(A19) 


in which the driving force is the concentration or density difference between the reaction 
zone electrolyte and the bulk electrolyte. Assuming that no diffusion takes place within 
the reaction zone, the average diffusive path length is one-half the thickness of the bulk 
zone (Zg/2). The densities in equation (A19) were calculated from a knowledge of the 
instantaneous zonal masses and their corresponding calculated volumes. The bulk 
volume and corresponding thickness Zg were calculated as outlined; the reaction zone 

volumes VV. and , which are assumed to be equal, were taken to be the instanta- 
H 2 °2 

neous volumetric sum of their respective reaction-zone agglomerate cylinders. 

The diffusion under consideration in equation (A19) takes place in a zone that is made 
up of both flooded electrode and matrix. For an assumed idealized pore structure, the 
development presented in reference 7 indicates that the tortuosity t is defined by the 
following function, 
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T _ 3(1 - a) 

n 


(M 


In this expression a. is the structure porosity and T / (3 the ratio of effective pore 
radius to solid fiber radius. By approximating the pore to solid radius ratio to be 2.0 
for both portions of the zone, tortuosity factors of 1. 05 for the electrode and 5. 0 for the 
matrix portions were obtained from reference 7. These data are presented as a family 
of linear characteristics that identify the tortuosity as being inversely proportional to 
structure porosity for various pore- radius to solid-radius ratios. Because of this differ- 
ence in tortuosities between the two sections of the bulk zone, the net tortuosity of the 
whole zone changes as the zonal volume varies. In the model this variable tortuosity is 
accounted for in the diffusion coefficient D R q-KOH - 111 P axticular > the individual tor- 

tuosities are treated as factors by which the liquid diffusivity D R o-KOH * s d * v i ded *° 

obtain the individual effective diffusivities for mass transport through the electrode and 
matrix portions of the zone. From these two values a weighted mean diffusivity 
D H o-KOH is determined for the overall zone; the weighting factor is determined from 

the instantaneous relative volumes of the electrode and matrix that constitute the bulk 
zone. 

Relating the zonal diffusivity to the time varying electrolyte volume in this manner 
represents the introduction of a system nonlinearity, which, along with the variable 
mass-transfer area, accounts for the experimentally observed, inlet humidity, nonlinear 
effect (fig. 6). 


THE GENERATION TERM 

The water generation rate w„ is represented in the basic conservation equation 

_ n sx _ n 

(eq. (Al)) by 1.87X10" (I) where the quantity 1. 87x10 kilograms per ampere per sec- 
ond is the Faraday constant. 

Therefore, as programmed for computer solution, the continuity equation for water 
in the hydrogen reaction zone is given by 


f 


M 


dM 1 


H 2°" H 2 

dt 


- 1.87X10 7 (I) - CA Z _ H (P - P g>i ) - d h o-KOH A Z 


H 2 °- H 2 


M 


h 2 °- b \ 


'Ho 


'B 


(A 20) 
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CONTINUITY EQUATION FOR WATER IN BULK ELECTROLYTE ZONE 

The general continuity equation for the bulk zone is 

dM, 


HgO-B 

dt 


- w e-h 2 ' w e-o 2 


where w E _ H represents the flow of water entering from the hydrogen reaction zone and 
2 

w E-0 flow exiting the bulk zone and entering the oxygen reaction zone. The ex- 
2 

pression for the flow from the hydrogen reaction zone into the bulk zone w E _ H was 

2 

developed for the hydrogen reaction zone conservation equation and is given by equation 
(A19). 

In applying Ficks' law to obtain w e _q , the driving force, which is the density dif- 

2 

ference between the bulk and oxygen zones, is the only parameter different from the ex- 
pression for w E _ E . Therefore, the expression for w e _q can be written directly as 


"HgO-B 


V. 


w E-0 2 “ d h 2 o-koh a z 


B 


m h 2 o-o 2 \ 

v o„ 


J B 


(A21) 


and the equation in programmed form is 




/ M H 2 °- H 2 

m h 2 o-b^ 


/ M h 2 o-b m h 2 o-o 2 \ 

dM H 9 0-B _ 

£ = Drr ^ ™ ti A„ 

' % 

v B 

_ Tj A 


dt h 2 o-koh Z 

1 

1 

n 2 a- koh a z 



V 2 

i 


V 2 / 


(A22) 
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CONTINUITY EQUATION FOR WATER IN THE OXYGEN REACTION ZONE 


The general conservation equation for the oxygen zone is given as 


dM 


Hij0“0« n 

— ± = W w r. - 0. 935X10" '(I) 

dt °2 


(A23) 


In this expression Wg_Q represents the flow of water entering the oxygen zone from 

2 _ rj 

the bulk zone, and 0.935x10" (I) the rate at which water is consumed by the oxygen- side 


reaction. The flow term w 


jg_Q is defined by equations (A21), and the numerical coef- 
2 ,-7\ 


ficient on the sink term (0. 935x10" ) is the Faraday constant. By applying the flooded 
agglomerate reaction zone concept and the ensuing assumptions as enumerated for the 
hydrogen reaction zone, the programmed form of the equation is given as 


/“ h 2 o-i 


M 


dM 


h 2 o-o 2 

dt 


" D H20-K0H A Z 


B 


H 2°-°2\ 


J B 


- 0.935X10" '(I) 


(A24) 


CONTINUITY EQUATION FOR OH" IN THE HYDROGEN REACTION ZONE 


As in the case of water movement, the transport of hydroxide ions dictated by the 
cell reactions can be mathematically simulated by the continuity equation. For the 
hydrogen zone it can be expressed as 

dM oir _ jj 

“ = W RZ ' 1 * 78x10 " 7 ( I ) ( A25 ) 

dt 

where w^g represents the transport of OH ions into the reaction zone and 1. 78x10 (I) 

the consmnption rate of OH" by the hydrogen reaction. The term defining the flow into 
the zone w^ i s analogous to the water flow terms outlined in the water continuity equa- 
tion except that the transport in this case results from ionic migration as well as diffu- 
sion. 

By assuming dilute solution theory to apply plus the absence of convective transport, 
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an effective mass-transfer coefficient can be approximated that accounts for both diffu- 
sion and ionic migration. The K + and OH" ionic fluxes are defined in a general sense 
by 


dC, 


N OH - " ( D OH") d ° H ' feir) ( U OH”) F ( 'OH') 


de 

dZ 


(A26) 



D 


dC 


K 


K + / dZ 



de 

dZ 


(A27) 


The first term to the right of the equal sign in each equation, Ficks' law, defines the 

ionic movement due to diffusion, and the second term defines the ionic movement due to 

electrical migration. Since the potassium does not enter into the electrode reactions, it 

can be assumed that the flux of potassium ions N is equal to zero; therefore, equa- 

K + 

tion (A27) can be solved for the potential gradient de/dZ in terms of the parameters 
that appear in equation (A26). By substituting this expression for potential gradient into 
equation (A26), the OH" flux equation becomes 


U, 


N, 


OH" 


= - ID, 


OH" 


dC 


OH 


U 


D 


OH 


K7 dZ 


r 


By the assumption of infinite dilution, which is inherent in the ideal solution approach, 
the ionic mobilities are defined by the Nernst- Einstein equation to be 


OH" ~ 


P OH" 

RT 


and 




Substituting these relations into the expression for Nqjj- determines the effective 
mass -transport coefficient D’ to be 2 Dqjj~. To a first approximation, the value 
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, binary diffusivity for the 


°f Uqjj may be represented by the value of Djj q_koh 
H gO-KOH system. ^ 

The transport equation for the movement of OH - in the hydrogen reaction zone can 
be expressed as 


A>H--B Moir - H 2\ 

dM OH--H, 

^ = D'A, 

dt 


CONTINUITY EQUATION FOR OH' IN THE BULK AND OXYGEN REACTION ZONES 

Except for the difference in parameter values and driving force densities, this equa- 
tion is analogous to the continuity equation for water in the bulk zone and can be written 
directly as 


V B 


Hr 


- 1. 78X10" '(I) 


(A28) 


dM, 


OH -B 
dt 


= D'A r 


/ lvl OH'-o 2 

m oh"-b\ 

v ° 2 

V B 




'M, 


M, 


OH“-B 


- D’A r 


OH - Hr 


Hr 


(A29) 


B 


\ -2 

V 2 


In this equation the mass of OH” contained in the oxygen reaction zone is determined, 
not from an independent continuity equation, but from a knowledge of the continuously 
calculated bulk and hydrogen zone OH - masses and the mass of pure KOH loaded into the 
cell at time of construction. Since the initial electrolyte loading was 0. 0105 kilogram 
and since a net removal of KOH is assumed not to occur, the mass of KOH present in the 
oxygen zone at any instant is 

m koh-o 2 = °' 0105 " m koh-h 2 " M KOH-B 


From the assumption of electroneutrality and the accompanying K + and OH - equality, 
the instantaneous oxygen zone OH" mass is determined from the preceding expression 
to be 
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+ 1 


- 0.0032 - M oh _ h ^ m qh"-B 


The largest error introduced by the electroneutrality assumption, as it is applied 
in this analysis, exists during periods of transient operation. However , the approxima- 
tion is justified during these periods by the fact that the ionic mobility is such that the 
realinement to neutrality is fast when compared to the transport transients. 


TRANSIENT TEMPERATURE EFFECTS 

In the type of cell under consideration, the effect that temperature has on the trans- 
port of water is manifest through the temperature- sensitive electrolyte vapor pressure. 
This sensitivity of vapor pressure to temperature is illustrated in figure 10 in which 



Figure 10. - Vapor pressure of aqueous potassium hydroxide as function of weight percent 
potassium hydroxide (ref. 8). 
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vapor pressure is plotted against weight- percent of KOH for a range of temperature 
values. 

In terms of water transport, a transient temperature increase that accompanies a 
load step produces an increase in electrolyte vapor pressure and hence an effective in- 
crease in the force driving water into the reactant stream. Therefore, unless compen- 
sated for by other cell operating parameters, an increase in cell temperature has a dry- 
ing effect on a cell. 

Since in the model the cell thermal response accompanying a step in load was not 
simulated analytically, experimental thermal transient data that were taken were used 
to calculate the time varying temperature effect any particular load transient had on the 
electrolyte vapor pressure. By knowing the approximate initial electrolyte concentration 
and cell temperature, a transient change in vapor pressure resulting from a transient 
temperature response was obtained by scaling the data of figure 10 for temperature and 
picking off data points to correspond to the data points of the time varying temperature. 
However, since the electrolyte weight fraction changes during any load transient, the 
vapor pressure data could not be taken off a constant weight fraction line in figure 10. 

To account for this fact in the model, an approximation for the final, posttransient con- 
centration was determined and the data points were picked off a straight line struck be- 
tween the initial and final concentrations. 

The transient vapor pressure increase resulting from the temperature transient 
accompanying a 40- ampere load step is shown in figure 11. The actual temperature cor- 



Figure 11. - Temperature effect on electrolyte vapor pressure for 40-ampere load step. Operating conditions: cur- 
rent step, 20 to 60 amperes,- inlet stream hydrogen content, 0.05 kilogram per hour (0. 12 lb/hr); water content, 
of inlet stream, 0.05 kilogram per hour (0. 12 lb/hr). 


39 



recti on was made by adding the time varying vapor pressure change of figure 11 to P* , 
the reaction zone electrolyte vapor pressure, as the model was responding the 
40- ampere step change. 


PERTINENT FUEL CELL PARAMETERS, DIMENSIONS, 
AND OPERATING CONDITIONS 


The following table lists values of pertinent fuel cell parameters, dimensions, and 


operating conditions. 

Measured values 

2 2 

Gas chamber (x direction) cross sectional area, cm ; in 5. 23; 2. 06 

Electrode area, A^, cm 2 ; in. 2 2 13 ; 33 . 1 

Hydrogen gas chamber length, L, cm; in 2. 21; 0.87 

Hydrogen supply chamber thickness, cm; in 0.152; 0.06 

Nominal cell operating temperature, T, K; °F 366; 200 

Compressed matrix thickness, cm; in 0.051; 0.02 

Catalyst loading, w, kg/cm 2 ; lb/in. 2 4. 0x10"^; 5. 68x10”^ 

Catalyst density, p, kg/cm 2 ; lb/in. 2 2. OlxlO” 2 ; 0. 724 

Electrode thickness, H, cm; in 0. 076; 0. 030 


Reasonable values obtained from the literature 

o 

Mass transport coefficient of OH" in electrolyte, D', cm /sec; 

in. 2 /sec 1 . 8x10” 4 ; 2. 8x10" 5 

Mass diffusivity of water - potassium hydroxide mixture, Djj q-KOH^ = D 

cm 2 /sec; in. 2 /sec 0. 90x10“^; 1. 4x10” ^ 

Bulk zone mean porosity, P 0.5 

Compressed matrix porosity 0. 42 


Assumed values 

Tortuosity of flooded electrode portion of bulk electrolyte zone 1. 05 

Tortuosity of matrix portion of bulk electrolyte zone 5.0 

- 4 

Agglomerate radius, r , cm; in 0.001; 3. 9x10 

Agglomerate porosity, e 0.50 
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APPENDIX B 


SYMBOLS 


A z 

a z-h 2 

c 

C OH" 

V 


D' 


^h 2 o-koh 

d h 2 o-koh 


D 


K* 


D, 


OH" 


2 2 

electrode area, cm ; in. 

2 2 

hydrogen reaction zone agglomerate surface area, cm ; in. 

2 

linear portion of overall mass-transfer coefficient, cm 

O 

molar concentration of hydroxide ions in KOH- water mixture, kg- mole/cm , 

O 

lb-mole/in. 

molar concentration of potassium ions in KOH- water mixture, kg-mole/ 

3 3 

cm , lb-mole/in. 

9 

mass- transport coefficient of hydroxide ions in cell electrolyte, cm /sec; 
in. 2 /sec 

2 2 

diffusion coefficient for H 2 0-H 2 , cm /sec, in. /sec 

O 9 

mass diffusivity of KOH- water mixture, cm'/sec; in. /sec 

2 

weighted mean diffusion coefficient for bulk transport zone, cm /sec; 
in. 2 /sec 

2 2 

diffusivity of potassium ions in KOH- water mixture, cm /sec; in. /sec 

2 2 

diffusivity of hydroxide ions in KOH- water mixture, cm /sec; in. /sec 


e electrostatic potential, V 

F Faraday's constant, C/equiv 

H electrode thickness, cm; in. 

2 2 

h local mass-transfer coefficient, mole/(cm )(sec)(N/cm ) 

hj local mass-transfer coefficient of transport interface controlled mole/ 

(cm 2 )(sec)(N/cm 2 ) 

h 2 local mass-transfer coefficient of transport diffusion controlled mole/ 

(cm 2 )(sec)(N/cm 2 ) 

I current, A 

2 2 

k. gas- liquid interfacial mass- transfer coefficient, mole/(cm )(sec)(N/cm ) 

L hydrogen gas chamber length, cm: in. 

I width of gas flow path in gas cavity, cm ; in. 
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m h 2 o-b 

m h 2 o-h 2 

M H2°-°2 

m koh-b 

M + 

K-B 

M OH"-H 2 

m oh”-b 

m oh"-o 2 

N 

•V 

n h 2 o 


p 

p 


h 2 o 


(Z) 


s,i 

K* 

^OH" 

R 


T 

t 


u 

U 

U 


OH“ 

K + 


mass of water contained in bulk electrolyte zone, kg: lbm 

mass of water contained in hydrogen reaction zone, kg; lbm 

mass of water contained in oxygen reaction zone, kg; lbm 

mass of pure potassium hydroxide in bulk electrolyte zone, kg: lbm 
mass of potassium ions in bulk electrolyte zone, kg; lbm 

mass of hydroxide ions in hydrogen reaction zone, kg; lbm 

mass of hydroxide ions in bulk zone kg; lbm 

mass of hydroxide ions in oxygen reaction zone, kg; lbm 

2 

number of agglomerate cylinders per cm of electrode 

flow rate of H 9 0 vapor out of electrode per unit electrode geometrical 
^ 2 

area, mole/(cm )(sec) 

9 

flow rate of HgO in Z-direction per unit gas flow area, mole/(sec)(cm ) 

flow rate of H 0 O into GFHR per unit gas- liquid interface area, mole/(sec) 
(cm 2 ) 

o 

local x-dependent electrolyte equilibrium vapor pressure, N/m ; psi 

2 

overall mean reaction zone electrolyte vapor pressure N/m ; psi 

9 

partial pressure of HgO at Z in GFHR, N/m ; psi 

2 2 

local HgO partial pressure in hydrogen supply chamber, N/m ; lb/in. 
partial pressure in hydrogen stream at cell inlet, N/m ; lb/in. 
valence of potassium ions 

valence of hydroxide ions 

water vapor gas constant, (cm)(lb)/(kg)(K); in. -lb/(lbm)(R) 
agglomerate radius, cm; in. 
temperature, K: R 
time, sec 

dimensionless parameter, hA z /(w^,/n) 

mobility of OH - ions, (cm ) (mole )/(J) (sec); (in. )(mole)/(in. -lb)(sec) 
mobility of K7 ions, (cm )(mole)/(J)(sec); (in. )(mole)/(in.-lb)(sec) 
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w E-H 2 

w h 2 ° 

w k 

w s 

w s 

Wrp 

X 

Z' 

Z B 

Z H 

e 

0 


9 

n 


p 

P E 


Ph 2° 

\f/(u) 

V 1 » ^ 2 ’ ^3 


3 3 

bulk electrolyte transport zone volume, cm ; in. 

3 3 

hydrogen reaction zone volume, cm ; in. 

3 3 

oxygen reaction zone volume, cm ; in. 

2 2 

catalyst loading, kg/cm ; lb/in. 

mass flow of water into bulk electrolyte zone, kg/sec; lb/sec 

molar flow rate of H 2 0 in gas cavity, mole/sec 

water generation rate, kg/sec; lb/sec 

mass flow rate of HgO into gas cavity, kg/sec; lb/sec 

molar flow rate of HgO into gas cavity, kg-mole/sec; lb-mole/sec 

total molar flow of gas through gas cavity, kg-mole/sec; lb-mole/sec 

distance from cell entrance in direction of flow, cm; in. 

distance (in gas phase) from agglomerate base, cm; in. 

bulk electrolyte transport zone thickness, cm; in. 

nonflooded portion of each electrode, cm; in. 

agglomerate microporosity 

gas- liquid interface area per unit electrode geometric volume, cm”*; in.”* 
gas filled cross-sectional area per unit geometrical area, dimensionless 

9 

total pressure, N/m ; psi 

3 3 

catalyst density, kg/cm ; lb/in. 

Q 

density of KOH- water mixture in bulk electrolyte transport zone, kg/cm ; 
lb/in. 3 

O O 

density of water in bulk electrolyte, kg/cm , lb/in. 
dimensionless overall mass- transfer coefficient 

alternate forms of , the dimensionless overall mass- transfer coefficient 
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